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ABSTRACT 

We present a new group field theory model, generalising the Boulatov model, which incorporates 
both 3-dimensional gravity and matter coupled to gravity. We show that the Feynman diagram 
amplitudes of this model are given by Riemannian quantum gravity spin foam amplitudes coupled 
to a scalar matter field. We briefly discuss the features of this model and its possible generalisations. 

o : 

O . !• INTRODUCTION 

<N ; 

Spin foam models Q, represent a purely combinatorial and algebraic implementation of the sum-over-histories 
approach to quantum gravity, in any signature and spacetime dimension, with an abstract 2-complex playing the role 
r , of a discrete spacetime, and algebraic data from the representation theory of the Lorentz group playing the role of 
geometric data assigned to it. Indeed, the first model of quantum gravity to be ever proposed, the Ponzano-Regge 
^vq \ model, was a spin foam model for Euclidean quantum gravity without cosmological constant Q. This approach has 
1 recently been developed to a great extent in the 3-dimensional case. It is now clear that it provides a full quantisation 
of pure gravity 0, whose relation with the one obtained by other approaches is well understood 0,0. Moreover, 
matter can be consistently included in the picture [J, Q , providing a link between spin foam models and effective field 
theory living on a non-commutative geometry. This picture allows us to naturally address the semi-classical limit 
of spin foam models and shows that quantum gravity in dimension 3 effectively follows the principle of the so-called 
deformed (or doubly) special relativity [9j . 

The group field theory formalism |lfj| represents a generalisation of matrix models of 2-dimensional quantum gravity 
■ • It is a universal structure lying behind any spin foam model for quantum gravity [l2l , providing a third 
quantisation point of view on gravity [T(j and allowing us to sum over pure quantum gravity amplitudes associated 
1 . with different topologies |14| . In this picture, spin foams, and thus spacetime itself, appear as (higher-dimensional 
^Jy analogues of) Feynman diagrams of a field theory defined on a group manifold and spin foam amplitudes are simply 
the Feynman amplitudes weighting the different graphs in the perturbative expansion of the quantum field theory. On 
the other hand, we can construct a non commutative field theory whose Feynman diagram amplitudes reproduce the 
coupling of matter fields to 3d quantum gravity for a trivial topology of spacetime |8j. Remarkably, the momenta of 
$J] ■ the fields are labelled also by group elements. Moreover, in three dimensions there is a duality between matter and 
geometry, and the insertion of matter can be understood as the insertion of a topological defect charged under the 
Poincare group 0. This suggests that one should be able to treat the third quantisation of gravity and the second 
quantisation of matter fields in one stroke (see 0] for an early attempt). The purpose of this paper is to study how 
the coupling of matter to quantum gravity is realised in the group field theory, and whether it is possible to write 
down a group field theory for gravity and particles that reproduces the amplitudes derived in coupling quantum 
matter to quantum geometry. 

This is what we achieve in the present work. The way the correct amplitudes are generated as Feynman amplitudes 
of the group field theory is highly non-trivial. It requires an extension of the usual group field theory (gft) formalism to 
a higher number of field variables, and produces an interesting intertwining of gravity and matter degrees of freedom, 
as we are going to discuss in the following. The formalism we present is still based on the classical SU(2), it bears 
strong similarity with a recent work of Krasnov 0] who considered gft based on the quantum group DSU(2). This 
should not be a surprise since it is well understood that the particle spin foam amplitudes have a quantum group 
structure hidden in them 0,0. The gft model we propose here is, however, very different from the ones considered in 
Krasnov's work since our Feynman graphs reproduce explicitly the insertion of particles coupled to gravity. 
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II. THE PONZANO-REGGE SPIN FOAM MODEL COUPLED TO POINT PARTICLES 



The general form of Feynman graph amplitudes for spinning particles coupled to 3 dimensional quantum gravity -eg 
the Ponzano-Regge model- has been written in Q . In this paper we focus on the case of spinless particles and we recall 
in this section the definition of these amplitudes before deriving them from a gft. We start from a triangulation A of 
our spacetime M and consider also the dual A*: dual vertices, edges and faces correspond respectively to tetrahedra, 
faces and edges of A. We choose our Feynman graph, "1, to be embedded in the triangulation A such that edges of "I 
are edges of the triangulation. Each edge of "I is labelled by an angle 9 E [0, tt] 

9 = nm, k — AttGn, 

where Gn is Newton's constant, n is the inverse Planck mass and m is the mass of the particle. To each angle 9 we 
associate an element of the Cartan subgroup H of SU(2) 

ie o 
e~ 1 ' 

which corresponds to a rotation of angle 29 around a given axis. Given a group G, here SU(2), we assign group 
elements g e * to all dual edges e* of the triangulation. We constrain the holonomies around dual faces /* ~ e to be 
flat if there is no particle and we constrain it to be in the conjugacy class 9 if e is an edge of "I. More precisely, let us 
denote by G e (— the product of the group elements around a dual face (or plaquette) /* ~ e: 



Ge = G f .= 1] 9. 

e*£df* 



where e/*(e*) = ±1 records the orientation of the edge e* in the boundary of the (dual) face /*. The amplitude is 
given by 



Z 



M 



O*) = A(0)I^I J l[dg e , [l[du e l[6(G e )l[5(G e u e heu- 1 ), (1) 



where dg is the normalised Haar measure and 5{g) the corresponding delta function on G, A(#) = sin(0) and \E-\\ is 
the number of edges in the particle graph 1. We see that two types of group elements arise in the construction of 
this amplitude, the g e * describe pure gravity excitations and the u e variables are associated with the particle degrees 
of freedom. They arise because the insertion of a particle locally breaks the Lorentz and translational symmetries of 
the gravity model and the former gauge transformation becomes dynamical at the location of the particle Q . The u e 
are then interpreted as giving the direction of the particle momenta propagating along the edge e. The fact that we 
are talking about spinless particles implies that u e should not be thought of as an element of G but as an element of 
G/H (with H = U(l)), that is momentum space. The insertion of spinning particles can be achieved by taking into 
account a non trivial dependence under the H part of u e . The main lesson which we learn from this amplitude, and 
which gives the key idea leading to a gft construction of such an amplitude, is the fact that we need both G variables 
describing the gravity excitation and G/H variables describing the propagation of particles. We can expand the 5 
functions in terms of characters 

Ha) = ^2d j Xj{g), 

3 

dj = 2j + 1 being the dimension of the spin j representation and perform the integration over <? e » , u e in order to 
obtain a state sum model 

{i.}eO eel t L 7et * Jet 5 Jet >> > 

where the summation is over all edges of A and the product of normalised 6j symbols is over all tetrahedra t. For each 
tetrahedron, the admissible triples of edges, e.g. (j eti ,j et ,j et3 ), corresponds to faces of this tetrahedra. Boulatov 
|l7j was the first to show that the amplitude JTJ can be obtained as a Feynman graph evaluation of a group field 
theory. It is important to note however that this amplitude is generically divergent. It is now understood IS] that 
this divergence is due to a translational gauge symmetry (equivalent on-shell to diffeomorphism symmetry) acting on 
the Ponzano-Regge model. This symmetry should be gauge-fixed in order to obtain well defined and triangulation 
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independent amplitudes. The gauge-fixing is easily implemented by choosing a maximal tree 1 T of A\~l. The gauge- 
fixed amplitude can be obtained from (yi by replacing the product J^[ e ^-| S(G e ) by a product over delta functions with 
e ^ T U 1. In terms of the state sum model Q the gauge-fixing inserts in the summation a factor YieeT^iefi which 
eliminates the sum over j e , e £ T. The overall amplitude does not depend on the choice of T. 



III. A GFT MODEL FOR 3D QUANTUM GRAVITY COUPLED TO SCALAR MATTER 

A. Action and Feynman rules 

We shall now define a field theory on a group manifold, whose Feynman expansion gives the above modified Ponzano- 
Regge model. We consider a generic real field <fi, over the Cartesian product of six copies of SU(2) 

52, 5 3 ; «i, «2, «8) : SU(2) x • • • x SU(2) -» M. (3) 

6 

This is the basic object of the theory and, just as in the other g roup field formulations of spin foam models, it has the 
interpretation of a '3rd quantised' chunk of quantum geometry |lCtll2|. However, in this extended formulation based 
on a 6-argument field, this chunk of quantum geometry carries also additional degrees of freedom, labelled by the 
extra u variables, that acquire the physical meaning of particle degrees of freedom (more precisely particle momenta) 
when a mass parameter is inserted in a suitable way, as we are going to show in the following. Let us now list the 
symmetries that this field is required to satisfy. 

• We require that 4> is invariant under (even) elements a of the permutation group of three elements S3, acting on 
pairs of field variables (<?i, Ui): 

<l>(9l,92, 33; Ui,U2,U 3 ) = 4>(g a {\) , g a (2), 9a(3)i Uo-(l), U<r(2)) <V(3))- ( 4 ) 



If we require the field to be invariant under even permutations of the three pairs of arguments, then this is 
equivalent to dealing with a complex field instead, with the odd permutations mapping the field to its complex 
conjugate |19j . and the Feynman amplitudes produced by the corresponding group field theory are in one-to-one 
correspondence with orientable 2-complexes, as explained in Il9| . We can more generally require the field 
to transform under an arbitrary representation (not necessarily reducible) of S3. This will affect the type of 
2-complexes generated by the perturbative expansion of the theory. We stress, however, that this would not 
imply any change for the amplitudes of the Feynman diagrams. 

• Furthermore, we pick a U(l) subgroup H, of SU(2), with the interpretation of the invariance subgroup for the 
particle momenta, and project three of the arguments into SU(2)/U(1) equivalence classes 

Pb$(gi,g2,g3;u 1 ,u 2 ,u 3 ) = T\dbi (f>(g 1 , 32,53; uih, u 2 b 2 , u 3 b 3 ), (5) 

so that the field becomes in fact a function over three copies each of SU(2) and SU(2)/U(1). 

• Finally, we project the first half of the field, i.e. the part dependent on the first three arguments, into its SU(2) 
invariant part, by imposing invariance under simultaneous right action of SU(2) on the first three arguments: 

P a (t>(gi,g2,g 3 ;ui,U2,u 3 ) = / da 4>(g\a, g 2 a, g 3 a;m, u 2 , u 3 ). (6) 

JSU(2) 

This last symmetry has a geometric interpretation, as in the usual Boulatov model. It imposes the closure of 
the triangle of which the field <p represents the 2nd quantisation, by constraining the spin variables dual to the 
gi variables associated to its three edges. 



A connected set of edges of A\~I passing through every vertex of A\~I 
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Given such a field, we can write down a Boulatov-like action, with the extra u variables simply mimicking the relations 
among the gravity degrees of freedom g: 



3 

dgidui[P a P b <f>(gi,g 2 ,g3; u lt u 2 , u 3 )] [P a Pb<l>(gi, 52, g 3 ; u x ,u 2 , u 3 )} 



A f 6 (7) 
+ 41 / II d 9idui [P a 1 Pb ± <j>(gi , g 2 , g 3 ; «i , u 2 , u 3 )] [P a2 P b2 0(54 , 55 , #3 ; u 4 , u 5 , "3 )] 



i=l 



x [P Q3 P &3 (j>(g 4 , 52 , 56 5 "4 , "2 , u 6 )] [P a4 P fc4 (/)(gfi , 55 , 5 6 ; mi , m 5 , u 6 )] . 



As we will see the Feynman amplitudes obtained from this model are proportional to those obtained by the Boulatov 
model, i.e. the usual Ponzano-Regge spin foam amplitudes describing pure 3d Riemannian quantum gravity. This 
shows that the u variables are completely redundant at this stage and do not have any real physical meaning. They 
are going to acquire it soon, however. Now we introduce a mass parameter in the theory, turning this redundant 
description of pure quantum gravity into a model for gravity coupled to scalar matter. We define a mass insertion 
operator Pg, acting on the field (j) as follows: 

Pe 0(51,52,53; U\, u 2 , u 3 ) = (t)(uihgu^ 1 gi,g 2 ,g 3 ;uie 1 u 2l u 3 ), (8) 

where hg — exp(9Jo) G H; 9 is half the deficit angle created by the presence of a mass m, 9 — 4nGm, j^=t being the 
Planck mass; J is the generator of the U(l) subgroup H, the same subgroup under which the Ui variables of the field 
are invariant and e is the non-trivial Weyl group element, given in the fundamental representation by: 



e 



-1 

1 



We define then a new group field theory model, representing 3d Riemannian quantum gravity coupled to scalar matter, 
whose dynamics are given by the action: 

S[4>] =^ J J\dgidui([P a Pb(l>(gi, 52, g 3 ;u 1 ,u 2 ,u 3 )][P a P- b (j)(g 1 , g 2 , g 3 ;ui,u 2 ,u 3 )] 



i=i 



A ' 6 



[P a P b (j)(gi, 52, 53; ui,u 2 , u 3 )][P e PaP- b (j){gi,g 2 , g 3 ; m, u 2 , u 3 )}) 



(9) 



4! 



F[ dgidui [P ai P bl <p(gi , g 2 , g 3 ; u x , u 2 , u 3 )] [P Q2 P b2 ^(54 ,g 5 ,g 3 ;u4,u 5 ,u 3 )} 



i=l 



x [Pa 3 Pb 3 0{g4,g 2 ,g6;u4 : ,u 2 ,u 6 )}[Pa 4 Pb 4 (j>(gi,g5,g6;ui,u 5 ,u 6 )], 



where a is a free parameter and all integrals are with respect to the normalised Haar measure. Let us write the 
Feynman rules of the theory in coordinate space. We have to identify kinetic and vertex operators. For this purpose, 
we write the action as 

3 3 

S[4>] =- / JJdgiduj Y[ dgjduj (j)(gi,Ui)IC(gi, gj,Ui,Uj)(j)(gj,Uj) 

i=i 3=1 ( 10 ) 

Y^dg^duij V(g ij , Uij )<t>{gij,u 2j )(t){g 2 j,u 2 j)4){g 3j ,u 3 j)4){g4 3 ,u4j). 

where in this integral, i ^ j, and (f)(gij,Uij) — 0(<7ii> g\ 2 , 513; tin, Ui 2 , Ui 3 ), and so forth. The kinetic and vertex 
operators are 

„ 3 „ 3 3 

fc{gi,9j,Ui,Uj) = / \\dbi 8(g l gl 1 )5{u i b i u^ 1 ) - a j J| db, L 8{u 1 h e u^ 1 g 1 g^ 1 )5{u 1 bieu^ 1 )^5{gjgj 1 )5{u j bjU^ 1 ), 
J i=i J i=l j= 2 

(11) 

\ f 4 

V(gij,Uij) = — / Y[ da, J[ dhj Siajg^gijOti ^{b^u^u^b^). (12) 
" J i=i j>i 
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We have purposefully discarded the a variables in the <5-functions of the kinetic term. This does not change, up to an 
overall factor, the computation of the amplitudes since Pg commutes with P a . We have also sidelined the sum over 
permutations, for ease of notation. Care should be taken in inverting the kinetic term, however, on the subspace of 
symmetric fields only. On this subspace, the kinetic term is indeed diagonal and we can proceed as follows. We now 
define the operators I and Kg: 

fC(gi,gj,Ui,Uj) = I - a Kg. (13) 

The propagator is the inverse of the kinetic operator. Furthermore, the operator Kg satisfies (Kg) 2 = I, as laid out 
below: 

(Kg) 2 = [ d 3 gd 3 u Kg(g,g,u,u)Kg(g,g,u,u) 

d 3 gd 3 ud 3 bd 3 b 5(ui^Mr 1 5i5r 1 )^(ff25^ 1 )^(ff3g3 _1 )^(wi^iewr 1 )^( w 2&2^ 1 )^(u3fe3'" 3 _1 ) ( 14 -* 



x 5(u 1 h g u 1 1 g 1 g 1 1 )S(g 2 g 2 ^(fib^ 1 )S(u 1 b 1 eu 1 l )8(u 2 b 2 u 2 1 )8(u 3 b 3 u 3 x ). 
We integrate with respect to the g variables 

(Kg) 2 = J d 3 ud 3 bd 3 bS(uihgu^\ 1 hgu^ 1 g 1 g^ 1 )5(g 2 g2 1 )S(g 3 g3 1 ) 

x 5(uibieu^ 1 )5(u2b2U2 1 )5(u 3 b3U^ l )5(uibieu^ 1 )5(u2^ 

and then the u variables 

(Kg) 2 = J d 3 bd 3 b5(u 1 b 1 ehee^ 1 b^ 1 u- 1 u 1 heu^ 1 g 1 g^ 1 )S(g2g 2 ~ 1 )S(g 3 g- 1 ) 

x 8(u 1 biebieu^ l )8(u2b2b2U2 1 )8(u 3 bzb 3 u^ 1 ). 



(15) 



(16) 



But b\ehge 1 b 1 1 = b\h e 1 b 1 1 = h g 1 since ehgt 1 = h g x and since b\ and h g 1 are in the same commutative U(l) 
subgroup. Furthermore eb\e = — b^ 1 . Finally, redefining (— bib^ 1 ) — > b\, 6262 - * &2 and 6363 — > 63 gives us 

(Kg) 2 = J d 3 b 8( 9l g~ ^fe^^feSs" 1 )^"!^^ 1 )^^^^ 1 )*^^^ 1 ) = /• (17) 

This leads to a nice closed form for the propagator 

V(gi,gj,Ui,v,j) = I ^ a J ■ ( 18 ) 



B. Feynman amplitudes and spin foam formulation 

To construct a generic Feynman amplitude, we will analyse the structure and gluing properties of the propagator 
and vertex operator: 

1. Vertex Operator 

We scrutinise (tT^ll in two parts: 

V(9ij,Uij) = — [ f[daiY[dbij 8(0 ,g,/ g,,n , : i.v/;,,,/,///,,//, i . 

4 - J 1=1 j>i " . ' 

g variable part u variable part 

The ^-functions over the g variables are the usual holonomies around the six wedges dual to the edges e, of a tetrahedron. 
So the model already has the structure of a 2-complex dual to a triangulation. The u variables represent the momenta 
of the particles and as such are identified with the edges of the tetrahedron. Each edge of the tetrahedron is shared 
by two triangles. The (5-functions above ensure that the momentum associated to an edge is the same when viewed 
from either of these triangles. This extra structure is not present in the Boulatov model. 
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2. Propagator 

The operator (|18l) glues two tetrahedra at a triangular interface. From the analysis of the vertex operator above, 
we know that each triangle of a tetrahedron has three wedges (dual to each of its three edges) and three momenta 
associated to it. The propagator has two terms with different gluing properties. 




For Vmassiess, the <5- functions over the g variables glue three wedges of one tetrahedron to three wedges of another 
tetrahedron pairwise, giving the holonomies around three composite wedges. The ^-functions over the u variables 
ensure that the momenta on the edges of the two triangles match. For V mass i ve , two of the 6- functions over the 
g- variables act as for V m assiess, similarly for the u variables. One of the <5-functions, however, couples the g\ and u\ 
variables, effectively placing a massive particle at the point where the edge of the tetrahedron intersects the composite 
wedge. The final ^-function to consider, 8{u\bxeui 1 ) 1 places a tag on the edge e, of the tetrahedron with the mass 
insertion. The tag ensures that if another propagator further around the sequence of Feynman graph edges e*, forming 
the boundary of a (dual) face /*, inserts a mass along the edge e, it will cancel with the first according to the property 
{Kg) 2 — I. We find in general that we can only have two possibilities when a face is fully assembled: no particle on 
the edge e, or one particle on the edge, according to whether there have been an even or an odd number of insertions 
of Kg on the boundary of the face. To be clearer, we will calculate more explicitly the amplitude for a generic dual 
face in the next subsection. In the end the partition function for the field theory, when expanded in Feynman graphs 
(i.e. in a power series for A) takes the form: 

Z= fvfie-W^-^Lzn (19) 
J ^ sym[T\ 

where v[T] is the number of vertices in the Feynman graph; sym[r] is its symmetry factor; and Z[T] is the amplitude 
for each Feynman graph, being given explicitly by: 

z\Y] = up] f n da e , n du e n s(q e ) n ^g^^u- 1 ), (20) 

where J\f[T] is a normalisation factor, discussed later, arising from the (partial) redundancy of the u variables in each 
diagram and the a dependence; a e * is the holonomy along an edge of the Feynman graph; G e = Tie*<zdf* a e^ 1 ^ s ^ ne 
holonomy around a face /* of the Feynman graph, which is dual to the edge of e of the triangulation; and "1 is the 
set of edges of the triangulation that have a particle present. We recognise in l(2l7|) the Ponzano-Regge amplitude Q). 
We see that the group field theory we have defined gives, in addition to a sum over all possible quantum gravity spin 
foams arising as usual as Feynman graphs of the theory, a sum over all possible massive spinless particle insertions 
in the spin foam, interpreted as a sum over all possible Feynman diagrams for a scalar field theory. Each gravity + 
particle configuration is weighted exactly by the amplitude of the Ponzano-Regge model coupled with massive spinless 
particles given in Q , and provides us also with a definite normalisation factor for each of these amplitudes. 



C. Amplitude for a generic face of the Feynman graph 

Consider a face /* of a Feynman graph, the boundary of which is formed by N contiguous edges, e*, labelled 
e*, . . . , e* N . The amplitude for this face is 

Mf*) = J(d... )vfvQ . . . r f ;v f N \, (21) 

where V( are the (^-functions from the propagator along e* relating to the face /*; VL are the (5-functions from the 
vertex where e* and e* meet, pertaining to the face /*. If we contract all the g variable (5-functions we get a final 
(5-function of the form 

N 

Mf*) = J]] / [d. .. )8(G f*u ai heu~^ x ■ ■ • x u ari hgu~^) x (5-functions over the u variables), (22) 

n=0 ^ 
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the sum is over different combinations of mass insertions; n < TV is the number of particle insertions in that specific 
term and Gt* is the holonomy around the face built up from a product of a e *. If there are n particle insertions in 
the g variable part then there will be n <5-functions in the u variable part with e inserted. Once we contract the u 
variables we get that the masses cancel just as in the calculation of (Kg) 2 = I. In the end, for each term in the sum 
we have two possibilities: For n even we get 

6(G f .)8(bf.), (23) 

where bf is a product of the b variables around the face. Therefore we get no particle insertion and a pure gravity 
face modulo an extra factor which we take into the normalisation. For n odd we get 

5(G/.ufcflU- 1 )<5(6/.e), (24) 
thus a single particle insertion and another factor which we take into the normalisation. 



D. Overall normalisation of the Feynman graph 



For a kinetic term with the structure I — aKg, the propagator takes the form: 

1 



1 - a 1 



(/ + aKg) , 



as we have shown. This produces a-dependent amplitudes when the expansion in Feynman graphs is performed. We 
have not specified what this a is and how it depends on the physical parameters of gravity or matter; indeed there is 
quite some freedom involved in choosing a specific expression for a and only further analysis of the model we proposed 
can narrow the range of possibilities down to a restricted one 2 . This parameter will enter the normalisation coefficients 
controlling the relative strength of Feynman diagrams with and without particles. The normalisation factor will clearly 

contain an overall factor (l — a 2 ) ' e ' where |e*| denotes the number of dual edges of the two complex. There will 
be an additional factor a each time Kg is inserted along a dual edge. If an even number of Kg are inserted along a 
face no particle circulates along that face. A useful formula in order to get the right normalisation factor in a given 
example is 

(/ + aK e ) n = i ((1 + a) n + (1 - a) n ) I + X - ((1 + a) n - (1 - a)™) Kg. (25) 

Along with this numerical factor, there is a singular factor coming from a redundant (^-function for each face as shown 
in (|23I24|) . For each face not carrying a particle we have a factor 

f dbS{b) = ^j + l), (26) 

Mi) 

and for each face carrying a particle we have a factor 

db5(be) = J2(2j + l)(-iy, (27) 

the sum, being over integer j, is obtained from the character expansion of the delta function. These expressions are 
unfortunately ill defined and they need a regularisation. A proper regularisation that can preserve all the symmetries 
of the theory is to replace the usual SU(2) group by a quantum group J7 g (SU(2)). It is expected that with this choice 
the key features of the model can be preserved and that the corresponding normalisation coefficients are given by 

±J2j + l],t* - {l _ q - H){l _ qH) , (28) 



2 a natural possibility in view of is to consider a = A (8) 
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with q = exp(^ 2 jv+i )' [ n ]"3 = (l n ~ Q~ n )/{<1 ~ I 1 ) an d £ = — 1 for a face with particle and t = +1 otherwise. Let 
us recall that if we consider the original Boulatov model such factors do not arise since they come from a redundant 
summation over spins dual to the u variables. The Feynman graph amplitudes of the Boulatov model are, however, 
not equal to the physical quantum gravity amplitudes. In order to get the quantum gravity amplitudes one has to 
divide out the infinite volume of a gauge symmetry. This is conveniently done by restricting the summation over spins 
as presented in the first section. This gauge-fixing procedure which is well defined at the level of the quantum gravity 
amplitudes is, however, not fully understood at the level of the gft. When we extend the gft to include the momenta 
variables u, we also extend in a trivial way the gauge symmetry of the quantum gravity amplitude. This is where 
the additional factors l|26l I27fl come from. The gauge-fixing is trivially realised by fixing to the spin dual to the 
variables u e . This is, we feel, one of the main open challenges in this domain: to understand if such a gauge fixing can 
be understood already at the gft level both for the Boulatov model and for our extension including particles; that is, 
whether we can already for the Boulatov model identify at the level of the gft the translational (or diffeomorphism) 
symmetry responsible for the divergences of the naive gravity amplitude. A similar identification should also be 
implemented for our particle model. We do not resolve this issue in the present work. 

IV. DISCUSSION 

A. Features of the model 

We have seen that the model correctly generates spin foam configurations with some dual faces carrying particle 
data, i.e. a mass label, indicating that a particle of the given mass is propagating along the edge of the triangulation 
dual to that face. This means that the group field theory produces all possible Feynman graphs for a scalar field 
embedded in the triangulation dual to the quantum gravity 2-complex, specifying the field propagator on each line of 
the Feynman graph, and this only. Interestingly, this is enough, in this 3-dimensional setting, for specifying fully the 
dynamics of the particles, i.e. their interaction. In fact, this is dictated by the Bianchi identity constraining the sum 
of curvatures in the boundary of any 3-cell of the dual complex around any given vertex of the triangulation. When 
one or more particles are meeting at that vertex, thus interacting there, this implies momentum conservation for their 
interaction, which is the only content of any <fi n theory. Because any number of particles can be incident to any given 
vertex of the triangulation in the model we propose, this means that this corresponds to a scalar field theory with a 
potential given by a sum over any power of the field operators: tp 3 (x) + (p 4 (x) + 

We have seen that the crucial property of the modified kinetic term we propose for producing mass insertions in the 
spin foam amplitudes is, besides the extension of the field to 6 arguments, the property (Kg) 2 = I for the operator 
Kg inserting the mass of the particles in the group field theory action. It is interesting to note that the presence of 
this operator which inserts particles also breaks a symmetry that the pure gravity model possesses. This bears some 
similarity with the fact that particles in 3d can be understood as defects breaking the translational symmetry of the 
theory without matter The symmetry is the following: Let's consider the transformation 

<fr{gi,92,g3;ui,u 2 ,u 3 ) -> 4>(v 1 g 1 ,g2,g3;wiUi,U2,u 3 ), (29) 

where vj., wi are arbitrary fixed group elements. This transformation is clearly a symmetry of the pure gravity action 
J7J|. This symmetry is, however, broken by the insertion of a mass term and only the transformation 

<t>(gi,g2,g3;ui,U2,u 3 ) -> (f>{vigi, g2,g3;viUi,U2,u 3 ), (30) 

preserves the action ©. 

B. A direct generalisation 

The model we presented in section IIIII produces mass insertions in different faces of the spin foam 2-complex by 
application of the operator Kg in the 1st argument of the field, carrying the gravity variable g\ . Because of permutation 
symmetry, the fact that one has chosen the 1st argument of the field for inserting a mass parameter and not, say, the 
2nd is irrelevant, as one can easily convince oneself. Still, one may find the fact that a mass parameter is inserted 
in only -one- of the arguments of the field a bit unsatisfactory, for symmetry reasons. Here we want to discuss 
briefly what happens when one relaxes this condition. The result is that one can write down a generalised version 
of the model presented above, that is, however, basically equivalent to it, and leads to the same type of graphs 
being generated. One can consider defining a generalised kinetic term with matter insertions, defined using a sum 
of operators Kg(l), Kg(2), Kg(3), each Kg(i) inserting a mass parameter in the i-th argument of field, thus having 
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as kinetic term an operator with the structure I — a(K$(l) + Kg{2) + Kg(3)). It is obvious that a model like this 
would generate exactly the same type of graphs and amplitudes as the one we have defined above. It is also easy to 
realise that, once such an operator is introduced, there will be graphs in which the insertion of a Kg(l) and a Kg(2), 
say, in different propagators would produce the same amplitude that would have been generated by the use in the 
kinetic term of an operator of the form K$(l, 2), i.e. an operator inserting a mass parameter in both the 1st and 2nd 
arguments of the field at once; and indeed one could generalise further the kinetic term to an operator of the form: 
I - a(K e (l) +K e (2) + Kg(3)) - b(K e (l,2) + K e (2, 3) + Kg(l, 3)). Carrying this line of reasoning even further, one is 
led to the kinetic term: 

K = I- a(K B (\) + K g {2) + K e (3)) - b(K e (l, 2) + K e (2, 3) + K g (l, 3)) - cK e (l, 2,3) = I- K totaU (31) 

where the Kg(l, 2, 3) is defined as the operator inserting a mass parameter in all the first three arguments of the field. 
Again, this generalised kinetic term leads to the same kind of Feynman graphs and amplitudes, as it is easy to verify. 
In fact the structure of the amplitudes is determined by the property Kg(i) 2 = 1, as we have explained, so it is enough 
for each mass insertion produced by the operators Kg(i,j) and Kg(i,j,k) to satisfy that property in order for the 
resulting amplitude to be of the form we have described. Of course, the normalisation factors for the amplitudes are 
going to be different from those of the simpler model presented in section HTD and it will depend in general on three 
different coupling constants a, b and c. This gives additional freedom that may well turn out to be useful in some 
situation. 

It is interesting to note that there are several choices of coupling constant that lead to further simplifications. For 
example, one can show that in order for the added terms to satisfy (Ktotai) 2 oc I, then one needs to choose 

= 0, 6 = 0, or b = 0, a + c = 0, (32) 

that is 

£ = 7-^(1,2,3), or K = I - a (Kg(l) + K e {2) + K e {3) - Kg(l, 2, 3)) . (33) 

Finally, a simple and highly symmetrical choice is 

K=(I- aKg(l))(I - aK e {2))(I - aK e (3)) 

= 1- a{K e (l) + K e (2) + K e (3)) + a 2 (K e (l, 2) + K e {2, 3) + K e (l, 3)) - a 3 K e (l, 2, 3). 

Since the (/ — aKg(i)) commute, we can easily compute the propagator 



(I + aK g (l))( I + aKg(2 ))(I + aKg{3)) 
(1-a 2 ) 3 



V = n ■ ( 35 ) 



It is not clear at the present stage, however, which specific properties one should ask the gft propagator to fulfill. 



C. Alternatives 



After having discussed some possible generalisations of the proposed model leading to very similar structures, we 
would like to discuss two genuine alternatives to it. One based on a much simpler action constructed inserting a 
mass parameter in the simplest way in the usual 3-argument field, leading however to a problematic structure for the 
Feynman amplitudes, and thus showing the need for the 6-argument extension on which we have based our model. 
The other keeps the same structure of the model presented in section ITTT1 but inserts the mass parameter by means of 
a modification of the interaction term in the gft action, instead of the kinetic one. As we will show, this modification is 
completely harmless. Consider first a three argument field, the same on which the Boulatov model is based, and insert 
the mass parameter hg and the velocities for the relevant particles Ui in the 1st argument of the field. The model we 
obtain is therefore realised by forgetting about the presence of the it.; variables in the extra slots of the generalised 
field used in the model presented in section ITTll The action is then: 

S[<f>] = \ J~\\ dgidui (^(3i:52, 53)0(ffi,52,53) - a (j)(g 1 ,g2,g3)<j){u 1 hgu^ 1 g 1 ,g2,g3)) 

(36) 

A f 

+ 41 I Y\. dgi 0(3i, 52, 53)0(54, 55, 53)</>(54, 52, 5e)0(5i, 55, 5e)- 



(=i 
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The Feynman rules for this action can be read out easily, and the amplitudes constructed in the usual way. The kinetic 
operator is: 

3 . 3 

fc{9i,9j) =Y[S(9i9i 1 ) - a J du x S(u 1 h e u 1 ~ 1 g 1 g 1 7 1 )Y[6(g l g~ 1 ) = I — aKg, (37) 

i=l i=2 

while the interaction operator is the usual Boulatov one. It is easy to see from the expression for the kinetic operator, 
that the mass-inserting operator does not satisfy any property like (Kg) 2 = I anymore. Instead, its square gives: 

r 3 

(Kg) 2 = / dut du 1 8(u 1 h e u 1 1 u 1 h e u 1 1 gig 1 x ) 8(g i g i 1 ). (38) 

•* i=2 

This makes the propagator much more complicated, being given by: 

V(g i ,g j )=I + Y J (*K e ) n , ( 39 ) 

n>0 

resulting in considerably more arduous Feynman diagrammatics. An example of an amplitude of this model, for a 
typical dual face is: 

A(f*) = 5{G f .u ai hgu-l X • • • x u an hgu~l), (40) 

for n less than the number of edges bounding the dual face /*. We see that there is no 'multiple mass cancellation' 
anymore, and we end up having multiple mass insertions on each face, in the typical Feynman diagram, each of which 
is associated to a different SU(2)/U(1) velocity element. There is a possible physical interpretation of the resulting 
configuration, which is in terms of multi-particle states. In other words, we would have more than one particle with a 
given mass associated to a dual face and so to an edge of the triangulation. We are not in the position of being able to 
exclude this interpretation, or definitely reject the model that generates these configurations altogether; however, we 
feel that such an interpretation is problematic for at least two reasons: 1) it would imply that more than one particle 
is propagating along the same link of the triangulation, so they would be located at the same point in the manifold 
and interact together with other (multi-)particles at the vertices of the triangulation; most important, 2) interpreting 
the particle configurations in the triangulation as Feynman graphs of some effective field theory would become much 
less straightforward than for the model proposed in section ITTT1 and applying any procedure to extract this effective 
field theory, like that used in [8j, to the amplitudes generated by this 'simplified' model would be quite cumbersome, if 
at all possible. We note that the same problem of 'multiple mass insertions' is generated by most other modifications 
of the structure of the model presented in section ITTT1 affecting the 3 extra arguments of the field, i.e. the tt, variables. 
Such modifications lead to losing the property (Kg) 2 = I which is responsible for 'mass cancellation' in the dual faces 
when the operator is inserted more than once, and that leads to the presence of only one mass in each dual face, 
and furthermore, to consistency with the interpretation of the mass-labelled graphs in the triangulation as Feynman 
graphs of a field theory. 

We conclude by mentioning instead a harmless modification of the model, that may even turn out to be useful in 
future studies. One can keep the structure of the field to be a function of 6 arguments, and one can keep the same 
form for the pure gravity action based on this 6 argument field, but choose to insert a mass parameter not in the 
kinetic term but in the interaction term of the group field theory action. The group field theory action would then 
be: 



9 / II dgi du% [ p c p b^>(gi ,92,93;u 1 ,u 2 ,u 3 )}[PaPi(j>(gi,g2,93;u 1 ,u2,u 3 )} 

1 J i=l 

X f 6 

+ 41 / II dg% du% [ Pa i Pb i ^ gi ' g2 ' g3 > Ul ' U2 ' m )1 [ p "2 Pb 2 <Kff4 , 35 , 93 ; Ui , u 5 , u 3 )} 



(41) 

x [ p a s Pb 3 ( t ) (9i,92,96;u A ,U2,u e )] [P ai P bi ^(g-L , g 5 , g 6 ; m , u 5 , u 6 )\ 

A f 6 

+ It / n d5idMl ^[ p eP ai Pb l (i)(9i,g2,g3\u 1 ,U2,u 3 )][P a2 P b2 (l>(g i ,g b ,g 3 \u i ,u bl u 3 )\ 
' "* i=i 

x [P a3 P b3 (f>(g4 , 32 , 56 ; "4 , u 2 , u e )} [P a4 P bi c/)(g 1 ,g 5 ,g 6 ;u 1 ,u 5 ,u e )}. 

As we have anticipated, however, it is straightforward to verify that this action leads to the same type of Feynman 
graphs and amplitudes that one gets instead by modifying the kinetic term, the only difference being that one would 
get a different normalisation factor, in front of each amplitude. 



11 



V. CONCLUSIONS 



We have defined a new group field theory for 3-dimcnsional Riemannian quantum gravity, and constructed its per- 
turbative expansion in Feynman diagrams. These diagrams correspond to spin foam 2-complexes describing quantum 
gravitational degrees of freedom dual to 3d triangulations, as in the Boulatov model, but they also carry additional 
labels which describe massive spinless particles propagating in the spacetime one reconstructs from the gravity degrees 
of freedom. They have the interpretation of Feynman graphs for a scalar field theory embedded in the triangulation 
representing spacetime. The amplitudes for these diagrams have exactly the form obtained in Q from classical con- 
siderations, i.e. are given by the Ponzano-Regge model coupled to massive spinless particles, and are shown in || to 
admit an effective (non-commutative) scalar field theory description, confirming the above physical interpretation. 

The model presented possesses some quite non-trivial and interesting features, that we highlighted in the paper, 
and that deserve further analysis. This represents a first step in a program of analysing the coupling of matter and 
gauge fields to quantum gravity in the group field theory approach. The next steps would be first of all the inclusion 
of spin degrees of freedom and the construction of a group field theory reproducing the amplitudes given in Q for 
massive spinning particles. Second, one should develop the study of the gft observables in the presence of particles 
and its relation with spin networks with open ends. Next, one would like to show how the effective field theory picture 
for the particle degrees of freedom can be obtained directly, and possibly in a simpler way from the group field theory 
formulation. Finally, the problem of the coupling of gauge fields and the description of their interaction with both 
gravity and matter fields, again in the group field theory formalism, should be tackled. 

To conclude, as we mention in the text, one of the most pressing issues for our model and the Boulatov model is 
to understand whether there are symmetries at the gft level that justify the gauge fixing needed at the level of the 
Feynman graphs to reproduce physical 3d gravity amplitudes. This is necessary in order to promote this type of gft 
to a fundamental model of three dimensional gravity coupled to matter. 
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